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Abstract 



S 

■ The anomalous scaling in the Kraichnan model of advection of the passive scalar 

by a random velocity field with non-smooth spatial behavior is traced down to 
the presence of slow resonance-type collective modes of the stochastic evolution of 
fluid trajectories. We show that the slow modes are organized into infinite multi- 
plets of descendants of the primary conserved modes. Their presence is linked to 
the non-deterministic behavior of the Lagrangian trajectories at high Reynolds 
numbers caused by the sensitive dependence on initial conditions within the vis- 
cous range where the velocity fields are more regular. Revisiting the Kraichnan 
model with smooth velocities we describe the explicit solution for the stationary 
state of the scalar. The properties of the probability distribution function of the 
smeared scalar in this state are related to a quantum mechanical problem involv- 
ing the Calogero-Sutherland Hamiltonian with a potential. 



1 Introduction 

One of the basic open problems in fully developed hydrodynamical turbulence is the 
understanding of the origin of observed violations of the Kolmogorov |T| scaling. The 
violations indicate presence of strong short-distance intermittency in the turbulent cas- 
cade, i.e. of frequent occurrence of large fluctuations on short distances. Recently some 
progress has been achieved in the understanding of the analogous problem for the pas- 
sive advection of a scalar quantity by a random velocity field. The scalar is known to 
exhibit strong short- distance intermittency even if such is absent in the velocity field. In 
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the simplest model of the passive scalar, due to Kraichnan ||, one assumes a Gaussian 
distribution of time-decorrelated and spatially non-smooth velocities. The anomalous 
scaling of the scalar in this model was related in references [§][§Jf| to zero modes of 
differential operators describing the stochastic evolution of the flow. In the present 
paper we elaborate on this idea showing that the short- distance intermittency of the 
scalar is due to the presence of slow collective modes in the otherwise super-diffusive 
evolution of the (quasi)-Lagrangian trajectories of fluid particles. We show that in the 
Kraichnan model the slow modes, reminiscent of resonances in multi-body problems, 
are organized into infinite multiplets of descendants with the zero modes playing the 
role of primary objects. This structure might indicate the presence of hidden infinite 
symmetries in the Kraichnan problem. 

The other important feature of the Lagrangian flow in non-smooth velocities is its 
intrinsically probabilistic character: the Lagrangian trajectories of the fluid particles 
behave randomly even in a fixed velocity field. This phenomenon appears to be closely 
related to the presence of the slow modes in the stochastic flow of fluid particles. In 
more realistic velocity fields which are regularized on the viscous scale the effective 
stochasticity of the fluid trajectories is due to their sensitive dependence on initial 
conditions on scales shorter than the viscous one. We expect both phenomena: the 
presence of resonant slow modes in the Lagrangian flow and the non-deterministic 
character of the fluid trajectories, to be present in more realistic high Reynolds number 
velocity ensembles and to be responsible for their intermittency. 

The version of the Kraichnan model with smooth velocity fields, relevant for the 
description of the distances smaller than the viscous scale, has been intensively studied, 
see and |7|] to jn|. We return to this case developing further the tools of harmonic 



analysis used first for this model in M and ||. These tools allow a fast calculation of 



the the Lyapunov exponents for the flow of fluid particles found first in [[J and [T(| . We 
also explicitly construct the stationary state of the scalar relating its functional Fourier 
transform to a certain Schrodinger operator on the symmetric space SL(d)/ SO(d) 
where d is the space dimension. In particular, we compute the exponential decay 
rate of the probability distribution functions p(8) of smeared scalar values obtaining in 
three dimensions (and above) a result different from that of || [|TIJ |L3| . The discrepancy 
is traced to the contribution of correlations of different (pairs of) fluid trajectories 
disregarded in M |TI| ||rj[ . For the exponential decay rate of the Fourier transform of 
p(9) our results reproduce fully the calculations of [PlHlOl and confirm their semiclassical 
interpretation proposed in [IT3 . 



The paper is organized as follows. In Sect. 2 we present the Kraichnan model and 
obtain its solution employing a path integral formalism. Sect. 3 recalls briefly the anal- 
ysis of || and establishing anomalous scaling of the scalar by perturbative analysis 
of the scaling zero modes of operators governing the flow. The physical interpretation 
of the zero modes as scaling structures conserved in mean is the subject of Sect. 4. Sect. 
5 discusses the collective slow modes of the random flow of fluid particles. The analytic 
origin of the slow modes is unraveled in more technical Sect. 6. Sect. 7 describes the 
intricacies of the probabilistic description of fluid trajectories. Finally, Sects. 8 and 9 
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study in detail the case of Kraichnan model with smooth velocity field elaborating on 
the earlier results of || and of [0] to |I3| . Appendix A explicitly analyzes the slow modes 
in the relative motion of two fluid particles in non-smooth velocity field. Appendix B 
contains some more details on the smooth velocity case related to the results of || and 

2 Kraichnan model of passive scalar 

Let us consider an advection of a scalar quantity T(t, x) (the temperature) in d space 
dimensions. The time evolution of T is governed by the linear equation 

d t T + v - VT - kAT = f (2.1) 

where v (t, x) is the incompressible (V • v = 0) velocity field of the advecting fluid, k is 
the diffusion constant and f(t,x) is a given source term. Denote by R(t,x;t ,x ) the 
solution of the homogeneous equation 

(d t + v • V - k A) R(t, t ) = (2.2) 

with the initial condition 

R(t ,x;t ,x ) = 6(x-x ). (2.3) 

We shall call R(t, x; to, Xq) the evolution kernel and the corresponding operator R(t, to) 
the evolution operator. The solution of Eq. (j2.1|) has the form 

T(t,x) = J R(t,x;t ,y)T(to,y)dy + jf J R(t,x; s,y) f(s,y) dsdy (2.4) 

with T(t Q ) being the initial configuration of T at time to- 

There exists a functional integral formula for the evolution kernel which, for suffi- 
ciently regular v, may be easily given a rigorous sense as an integral with respect to 
the Wiener measure with density: 

t-./. . \ f —-r- f* ds \x(s)—v(s,x(s))] 2 „ /r _ _ N 

R{t,x;t Q ,x ) = / e iK Jt o lv ' v v ,n Dx (2.5) 

a;(t) — x 

where x = ^| . It will be useful to rewrite the above functional integral as a phase space 
one: 

D , \ f - f. ds[K P (s) 2 + ip{s)-(x(s)-v(s,x(s)))] , v 

R(t,x;to,Xo) = e Jt o Vx Dp (2.6) 

x(tQ) = XQ 

x(t) — x 

with the Gaussian integral over the unconstrained paths s i— > p(s) reproducing the 
previous integral. From the functional integral representations it is clear that when 
k — > then 

R(t, x; t , x Q ) 5(x - x(t;t ,x )) (2.7) 
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where x(t; t , xq) is the Lagrangian trajectory of the fluid particle satisfying the equa- 
tions 

X = v(t,x), x(to) = Xq. (2.8) 

Indeed, we may set k = in the phase space integral (|2.6| ) and the p-integral gives 
then a delta function(al) concentrated on the Lagrangian trajectory. For small positive 
k, on the other hand, R(t,x;t ,x ) is the probability distribution function (p.d.f.) 
of the endpoint of a small Brownian motion around the Lagrangian trajectory. Such 
a Brownian motion Xp(t\ t , Xq) with a drift is a solution of the stochastic ordinary 
differential equation (ODE) 

dx = v(t,x)dt + nd/3 , x(t ) = x (2.9) 

with f3(t) denoting the Brownian motion without drift. Thus 

R(t,x;t ,x ) = E(S(x - X/3(t;to,x ))) (2.10) 

where E( ■ ) denotes the expectation with respect to the Wiener measure of (3. Eq. 
(|2.10|) is another form of Eq. ( |2.5|) . 

We shall be interested in the situation when both velocities v and source / in 
Eq. ( |2.1| ) are random so that Eq. (|2.1| ) is a stochastic PDE. In order to solve such a 
stochastic equation, we should define the evolution kernel R(t,x;t ,xo) as a random, 
v-dependent process or, in plain English, be able to compute various expectation values 
of R'b like 

n n 

P n (t,x;t ,x ) = (JJi2(* i ,x i ;* ,i,a;o,i)) = E { II <K X * ~ x ^i t ^ *o,i, Xo,i)) (2.11) 

i=0 i=0 

where t = (ti, ...,t n ), x = (xi, . . . ,x n ) and similarly for t , x and where E(-) 
denotes the expectation w.r.t. /3j's and v. It is clear from the second expression for 
V n (t, x; t , Xq) that it gives the joint p.d.f. of the ends xp.(ti', to,j> Xo,i) of n Brownian 
motions (independent for given v) around the Lagrangian trajectories starting at times 
to t i from points Xq^. The k — > limit of V n (t, x; t , x ) (if exists) should simply give 
the joint p.d.f. of the endpoints x(ti] to,i, Xq^) of n Lagrangian trajectories 1 . 

Let us assume that the velocity is a Gaussian stationary field with mean zero and 
covariance 

(v a (t 1 ,x 1 )v p (t 2 ,x 2 )) = D a P(t 12 ,x 12 ) (2.12) 

where t i2 = ii — t%, x i2 = X\ — x 2 and d a D al3 (t,x) = to assure the incompressibil- 
ity Employing the phase space path integral representation (|2.6| ) and performing the 
Gaussian functional integration over v, we obtainF]: 

-T, I. . \ if - J2dftl dSi i K Pi( S i) 2 +iPi(Si)-(±i(Si)-v(s i ,X i (s i )))] } . 

V n [t, x; t , Xo) = ( e ^ Jt w DxDp ) 

x i i.^i) = x i 

1 K.G. thanks Ya. Sinai for attracting his attention to the statistics of Lagrangian trajectories 
2 similar expressions appeared in [|l4| |H| 
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I t *. dSi l K Pi( s i) 2 + i Pi( s i)'*i( s i)] } 



®i(*0,i) =!C 0,» 

31 £ 31 i 



e 2Zj «J% lJ %- 3 Vl J ' 8V w 3 3 1 J 3 DxDp. (2.13) 



If, following Kraichnan M } we assume that v(t, x) is also decorrelated in time, i.e. that 

D af3 {t,x) = 5(t)V a/3 (x), (2.14) 

then formula ( |2. 13| ) for V n (t, x; to, x ) may be further simplified. Let as set all U equal 
to t and all t 0i equal to t . We shall denote the corresponding V n by V n (t, x; t , Xq) 
(the general case can be reconstructed from the special one for u's delta-correlated in 
time). The fundamental property of the p.d.f.'s V n (t, x; i , x o) for the time-decorrelated 
velocities (not necessarily Gaussian) is the composition property 

Jv n (t,x; s,y)P n (s,y;to,x ) dy = V n {t, x; t , x ) . (2.15) 

From expression ( |2.13| ) we obtain, assuming relation ( 2.14|) , 
V n (t, x; to, x ) 

= /•g-^/^E^C^ + IE*,^^^ DxDp. (2.16) 

x £( £ 0,i)^ x 0,£ 
(*£)~ x £ 

It is easy to see that the right hand side is a phase space path integral expression for the 
heat kernel (dynamical Green function) of the 2 nd order (positive, elliptic) differential 
operator 



M„ = 4E V^(x i:j )d x? d x p -^A,, (2.17) 

i,j=l ' i=l 



i.e. that 

P n (t,x;t ,Xo) = e-^^^x^o), (2.18) 



compare Eq. (|2.6| ). Note that, due to incompressibility, there is no ordering ambiguity 
in passing from the path integral to the expression for M. n - Rigorously minded person 
may take expressions ( [2.18] ) as defining the evolution operators for the stochastic PDE 
equation (|2~T| ) in the time decorrelated case. Of course, the composition property (|2.15| ) 
follows from the semigroup law for the heat kernels. 

Let us now go back to the passive scalar. Assume that both v and / are independent 
stationary processes. Imposing also the zero initial condition for T at to = _ oo, we 
obtain using Eq. ( |2.4|) the following expression for the correlators of T: 



n n .tj . n 

Y[T{U, Xi )) = 1] / dSi d Vi P n (t,x;s,y) tfl ( 2 -!9) 

i=0 i=l J - aa -- 1 
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It should be clear that if a stationary state of the scalar is generated for large time and 
independent of the initial conditions (say, decaying at infinity) then its correlation func- 
tions should be given by the above equation. Hence the importance of understanding 
the behavior of the p.d.f.'s V n (t, x; t , x ) . 

Assume now that the source / (independent of v ) is a Gaussian process with mean 
zero and covariance 

(f(t 1 ,x 1 )f(t 2 ,x 2 )) = 6{t l2 )C{x 12 ) (2.20) 

where C is a positive definite test function. In this case and for the Gaussian, time 
decorrelated velocities, Eqs. ( |2.19| ) simplify permitting an inductive calculation of the 
stationary correlation functions of the scalar. Let us see how this works for equal 
time correlators. We may consider only the even-point functions of T, jF 2n (x) = 
(T(t,xi) ■ ■ ■ T(t,x 2n )), since the odd correlators of / vanish implying the same prop- 
erty of the T correlators. For the 2-point function we obtain 

T 2 {x 12 ) = [ dsfdy e-^- M2 (x,y)C(y 12 ) = / dy M 2 \x, y) C(y 12 ) (2.21) 

and for the 4-point function 

t s 

JF 4 (x) = ]T / dsf ds' fdy [dz e-^ M *(x,y)C(y t ,) 

•e" MA<, (l(i,^W,«)CN 

i 3 

= E / Mz\x,y) F 2 (yi,.^,yd C( Vij ) dy. (2.22) 

l<i<j<4 % 3 

Similar arguments for the 2n-point function give 

^2n(x) = M 2 ^(x,y)J 7 2n _ 2 (y 1 ,.^.,y 2n )C(y i:j )dy. (2.23) 

l<«<j<2n i i 

The above equations permit an inductive calculation of the stationary equal time cor- 
relation functions of T with the use of the (static) Green functions A^~ 1 (x, y) of 
operators A4 n . 

3 Zero mode dominance 

We shall be interested in the case where the spatial part T> a/3 (x) of the f-covariance 
has the form 

V aP {x) = V 5 a ? - d a ?{x) (3.1) 
where d al3 (x) scales with power 2 — 7, 

d aP (x) ~ |x| 2 " 7 , (3.2) 
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for small \x\. Here < 7 < 2 is a fixed parameter. The tensorial form of d a/3 (x) is 
fixed for small \x\ by the incompressibility condition d a d a ^(x) = 0: 

d a P(x) S l-7)5 a/3 |x| 2 - 7 - {2- 1 )x a x^\x\^) = dg(x) (3.3) 



where D is a constant. For < 7 < 2, one may take 



where m is an infrared regulator. Relations (3.2) and (|3.3| ) hold then for m\x\ <C 1. 
When m — > 0, c? a/3 (x) tends to the scaling form df@(x) but V diverges like 0( 



-ik-x 



k a k^ 



w?- 2 ). 

The Gaussian distribution with covariance given by Eqs. ( |2.14j ) and fl3.1| ) is relatively 
far from a realistic description of the statistics of turbulent flows. First, it excludes 
the velocity intermittency, i.e. more frequent occurrence of large deviations of velocity 
differences than in the normal distribution. Such occurrence characterizes short scales 
in the inertial interval of the turbulent cascade. Second, the time decorrelation is a 
brutal approximation since one observes scale-dependent time correlations in turbulent 
flows. The power-law growth of the velocity difference covariance 

((v a (t 1 ,x 1 )-v a (t 1 ,x 2 ))(vP(t 2 ,x 1 )- V P(h,x 2 ))) = 28(t 12 )d^(x 12 ) 

~ Sit^lxul 2 - 1 (3.5) 

mimics, however, the expected behavior in the turbulent cascade (the Kolmogorov value 
of the scaling exponent corresponds to 7 = | since time appears to scale like length 
to power 7 in the model). The point is that even the velocity distributions far from 
realistic, as the one described above, induce strongly intermittent scalar distributions 
and the purpose of the Kraichnan model is to understand this phenomenon in the 
simplest context. 

Operators M. n may be rewritten in the form 

n n 

M n = J2 d aP (xij)d X fd x g - kJ^A^ - jPoE^) 3 (3.6) 

l<i<j<n i=l i=l 

where the last term drops out in the action on translationally invariant functions. 
We shall, somewhat pedantically, denote the operator Ai n acting in the translational 
invariant sector by M n . We shall view M n as an operator in the reduced space L 2 (R dn ) , 
with d n = (n — l)d. This is the space of functions of the difference variables Xi n = 
Xi — x n , square-integrable with the measure d'x. = dx\ n - ■ -dx( n -i) n - The heat kernel 
of M, 



n ■ 



e -(t-to)M n / 



x ) = / "P n (t,x + a;t ,x ) da = P n (t, x; t , Xo) , (3.7) 



with a = (a, . . . , a) , gives the joint p.d.f. of the differences Xi n of the Lagrangian 
trajectories starting at points x (or, equivalently, the joint p.d.f. of the Lagrangian 



7 



trajectories in the quasi-Lagrangian picture jl6fl). It is translationally invariant sepa- 
rately in x and x . 

In the limit m — > when d al3 (x) takes the scaling form ( |3.3| ) but T>o diverges, 
operator M n , unlike A4 n , tends to the limit which for k = coincides with the scaling 
operator 

M s : = £ Mr* (3.8) 

\<i<j<n 3 

of scaling dimension —7. M^ c is a positive singular elliptic differential operator of the 
2 nd order in L 2 {H dn ). By a simple self-consistent analysis one may convince oneself that, 
at least for 7 close to 2, e~* Mn (x, y) and M n " 1 (x, y) converge pointwise when m —* 
and /t — > to the heat kernel e~* M ™ c (x, y) and the Green function (M^ c ) _1 (x, y) , 
respectively. The latter should satisfy bounds that may be inferred from a semi-classical 
analysis of the path integral expressions (p.!6|) with T> al3 replaced by rf^f and k set 



to zero^. In the limit 7—^2, d^(x) tends for non-zero x to a constant times 5 al3 
and M® c becomes proportional to the d n -dimensional Laplacian. When 7 is close to 
2, the heat kernel e~* M ™ c (x, y) and the Green function (M® c ) _1 (x,y) differ little from 
the heat kernel and the Green function of the Laplacian. In particular, e _tM ™ c (x, y) is 
finite everywhere and (M^ c ) _1 (x, y) is infinite only when x = y. We expect this to 
hold for all 7 > 0. When 2 — 7 is small, the behaviors of (M^ c ) _1 (x, y) around x = y 
and at infinity differ from those of the Green function of the <i ra -dimensional Laplacian 
by 0(2 — 7) modifications of the power laws. All that implies that the pointwise limits 
m — > and k —>■ of the equal time correlators of T given by Eqs. (|2.23j) exist, at 



least for 7 close to 2, and are given by the version of the same equations employing 
the scaling Green functions (M^ c ) _1 (x, y) (with dy replaced by d'y). From now on, 
we shall deal only with these limits and with the scaling operators M^ c and shall drop 
the superscript " sc" . 

We are interested in the behavior of the equal time correlators of T, especially in 
their scaling properties, in the situation when the source acts only on large distances, 
i.e., in our Gaussian model, when the spatial part C of the covariance of / is almost 
constant. We may study this regime by replacing C(x) by Cl(x) = C(x/L) and by ex- 
amining the large L behavior of the equal time correlators T^ n = J-2 n ,L ■ The following 
result, which may be referred to as the zero mode dominance, has been described in 
and : at 7 sufficiently close to 2 and at m, k = , 

^(x) = AcIfJ^W + 0{L- 2+0 ^) + [...] (3.9) 

for n > 1 . Above, Ac is a non-universal amplitude (a constant depending on the shape 
of the source covariance C) and p2n = ^pr^(2 — 7) + 0((2 — 7) 2 ) is a universal (i.e. 
C-independent) anomalous exponent. T^ n is the scaling (translationally invariant) zero 
mode of M 2n , 

J=l n (M) = A^-'-^Jx), M 2n Tl n = 0. (3.10) 



3 to our knowledge, such bounds have not been obtained in the mathematical literature and consti- 
tute an open mathematical problem 



s 



[ . . . ] denotes terms which do not depend on at least one Xi and as such do not 
contribute to the correlation functions of scalar differences (Yl(T(t,Xi) — T(t,yi))). 

i 

J- 2 ° n (x) = S x 2 12 x 2 u ■ ■ ■ xti, 2n + 0(2-7) + [ • • • ] (3.H) 

where S is the symmetrization operator. The contribution to jF° n proportional to 
2 — 7 is also known up to [ ■ ■ ■ ] terms ||. A similar analysis was performed in |4j] and 
[P7| for large space dimensions d. 



The main implication of the relation (|3.9|) is the anomalous scaling of the n > 1 , 7 
close to 2 structure functions S2u,l{x) = ((T(t,x) — T(t,0)) 2n ) . At m, k = and for 
\x\/L < 1, 

S 2 u,l{ X ) ~ L^|x| 7n - p2 ". (3.12) 

The above behavior contradicts the simple dimensional prediction 5*2 n (x) ~ |x| 7n 
which holds only for the 2-point function. 



Let us sketch the argument leading to the result ( [3.9D , based on applying the Mellin 
transform to select the dominant contributions for large L. It will be convenient to 
work with a version of operators M n of scaling dimension zero 

N n = Wj 2 M n Wj 2 (3.13) 

where R 2 = J2 t<3 ( x i ~ Xj) 2 - N n is also a positive (unbounded) operator^ in L 2 (R dn ). 
Since it commutes with the self-adjoint generator of dilations 

i 

it is partially diagonalized by the Mellin transform of the translationally invariant func- 
tions 

00 

/(x) ^ /(a,x) = [x-^fiXS^dX. (3.15) 



The map ( 3.15|) is a unitary transformation, diagonalizing D n , between 
L 2 (R dn ) and L 2 ({Rea = — ^-}) <8> L 2 (S dn ~ l ) , 

where S dn ~ l is composed of points x = x/R n in the space R dn of difference variables. 
In the language of the Mellin transform, N n becomes a family N n (a) of operators in 
L 2 (S dn ~ l ). In particular, 

(iV-V)>,x) = / N^(a;%y)f(a,y)dy. (3.16) 

Sdn-l 



4 technically, N n , as well as M n , may be defined as the Friedrichs extension of its restriction to 
smooth functions with compact support, vanishing around the diagonals Xi = Xj 
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where the Mellin-transformed Green function N n 1 (cr; x, y) satisfies the hermiticity re- 
lation 



iV-V;y,x) = N-\-d n -a;%y). (3.17) 
It is a meromorphic function of a with simple poles for generic 7. Around the poles 



N~\a- l;x,y) - _L_ /,(*) ft (y) (3.18) 

where /j are the scaling zero modes of M n of scaling dimension Oi and g>j are similar 
modes with scaling dimensions — d n + 7 — cfj, both in L 2 (S dn ~ 1 ) . Although operator 
M n has continuous spectrum when considered as a positive operator in L 2 (R dn ), it 
induces an operator N n (a — }) in L 2 (S' dn_1 ) with a discrete spectrum when acting on 
scaling functions with a scaling dimension a. The scaling zero modes occur at discrete 
values (jj of a for which zero belongs to the spectrum. 



It is easy to see from the inductive equations ( 2.23 ) that 

T 2n>L {yi) = L"^ 2n>1 (x/L) (3.19) 
and that, with the use of the Mellin transform, these equations may be rewritten as 

^,l(x) = L^j^r {Rn/LY N-\a-l x,y) 

• {^2n-2,x®CT{a -7, y) dy. (3.20) 
Shifting the integration contour to Recr = + 2 — 0(2 — 7) , we obtain 

^,l(x) = -E^^^/Res^iV-^a-Jx.y) 

■ (^-2,1 ® C)"(a - 7, y) rfy + 0(L- 2+ °( 2 ^)) (3.21) 

where the sum runs over the poles Oi in the strip 

-^t + l < Re^ < 7 n + 2 - C(2 - 7) (3.22) 

and the last term, suppressed for large L, comes from the shifted contour. There are 
two types of poles: those coming from (J r 2n-2,i ® Cf and those in the Green function 
N" 1 . The first ones contribute either to [ . . . ] or to £((£-2+0(2-7)) j n Eq. (gj) and 
are not interesting for us (at least for 7 close to 2). The second ones are related to the 
scaling zero modes of M n , see Eq. ( [3.18 ). Only rotationally invariant (if C has the same 
property) zero modes symmetric under permutations of points and square-integrable 
on S"^ 1-1 contribute to Ti n ,L- Such zero modes may be studied for 7 close to 2 by 
perturbative analysis of discrete-spectrum operators M n acting on scaling functions 
or, equivalently, of operators N n (a) acting in L 2 (S dn ^ 1 ). (Recall that for 7 = 2, M n 
becomes the <i n -dimensional Laplacian). The result is that, for 7 close to 2, all but 
one zero modes in the strip (|3.22|) contribute [ . . . ] terms. T 2n is the exception and 
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it has scaling dimension a = — p 2ri . We expect essentially the same picture with 
the zero mode domination of correlation functions to persist for all 7 > 0. One of the 
open problems is whether there are other non-[ . . . ] zero modes entering the strip 
( |3.22| ) for smaller 7 and whether, if they cross, they may produce pairs of zero modes 
with complex scaling dimensions. For 7 = the singularities in the inverse symbols 
of operators M n become strong enough to induce continuous spectrum of operators 
N n (a) and the picture of zero mode dominance has to be somewhat modified @||12] . 



One may also read the 0(2 — 7) contribution to the anomalous exponent p2 n from 
the 0((2 — 7) InL) term in the expansion of J^n.L into powers of 2 — 7, similarly 
as in the e-expansion for critical phenomena one obtains anomalous exponents from 
logarithmic divergences. In the latter case, the renormalization group which exponen- 
tiates the divergent logarithms provides an explanation why it is reasonable to extract 
information from badly divergent expansions. In our argument, the Mellin transform 
analysis played a similar role exponentiating the logarithms of L . One may show 



that there is an (inverse) renormalization group picture of the advection problem hidden 
behind the above analogy. The renormalization group for the passive scalar eliminates 
inductively the long-distance modes, unlike in critical phenomena where it is based on 
subsequent elimination of the short- distance degrees of freedom. 



4 Conserved scaling structures 

In view of the domination of the equal time correlators of the scalar by the scaling zero 
modes of operators M n , it is important to understand the physical interpretation of 
such modes. It is, in fact, very simple: 

zero modes are scaling structures preserved in mean by the flow. 

Indeed. Recall that e~ tMn (x, xo) = P ra (t,x;0,xo) and it describes the probability that 
the differences of n Lagrangian trajectories starting at time from points xo are at 
time t equal to Xi n . The mean value of a translationally invariant function /(x) of 
positions of n fluid particles at time t is then equal to 

(/>«,„ = / /(x)P B (f,x;0,x o )d'x = //(x)e- tM "(x,x Kx. (4.1) 

Differentiating the right hand side w.r.t. t, we obtain 

J /(x) M n e-* M "(x,x ) d'x = J M n /(x) e-* M "(x,x ) d'x (4.2) 

where we have integrated by parts twice. If / is a zero mode of M n then the right 
hand side vanishes and, consequently, the mean ( |4.1| ) is constant and 

(/U, = /(xo)- (4.3) 

In fact, the story is a little bit more complicated. The zero modes with Re<7j > 
— d n + 7 are true zero modes. However the ones with the real part of their dimension 
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< —d n + 7 are not. For them, M n f is a contact term supported at the origin. Such 
contact terms may give non-zero contributions to the right hand side of Eq. (|4.2|) or 
to the boundary terms in the integration by parts, depending on the interpretation. 
The zero modes with the scaling dimensions belonging to the strip ( |3.22| ) are true zero 
modes and hence they describe scaling structures of the flow conserved in mean. As was 
mentioned before, the translationally invariant zero modes that are square-integrable 
on S^ 1 come in pairs (fi,gi) corresponding to scaling dimensions cr, and — d n + ^ — Wi 
(we may assume that Re<jj > — + ^). For 7 close to 2 there are no zero modes 
square integrable on S^™ -1 with dimensions in the strip —d n + 7 < Re a < 0. We 
expect this to hold for any 7 > 0. In that situation are the true zero modes and 
they have non-negative real parts of the scaling dimension whereas are the false zero 
modes with real parts of dimension < — d n + 7 and with M n gi being contact terms. 
Our claim about the absence of zero modes in the strip —d n + 7 < Re a < may seem 
paradoxical if we recall that the multi-body structure of operators M n assures that 
zero modes of M n _ p are also annihilated by M n . IndeedQ, the (false) zero modes of 
M n _ p with Recr < —d n - p + 7 may lie in the strip — + J < Re<r < 0. However, the 
resulting zero modes of M n are not in L 2 (S dn ~ l ) and do not contribute to the poles of 
the Green function A^"~ 1 (x, y) and hence to the right hand side of Eq. ( |3.21| ). 



It should be stressed that the behavior (|4.3|) is atypical. For a general translationally 
invariant, scaling function with (say, positive) dimension a and for any time r > 0, 

(/>*,, = & J /(*) e"™"(x, (j) 4 x ) d'x (4.4) 
as it is easy to see with the use of the scaling property 

e- A7 * M "(Ax,Ax ) = A- d "e-* M "(x )X o). (4.5) 
It follows that, typically, 

</>, Xf) ~ t\ (4.6) 

The behavior (03) characterizes a super-diffusion where the square distances be- 



tween points grow faster than linearly in time. A slower behavior requires vanishing of 
J /(x) e _rMn (x, 0) (fx. Note that the exponent of the growth diverges when 7 — ► 0. 

It is easy to understand the origin of the behavior ( |4.6| ). The stochastic process 
described by the probabilities P n (t, x; 0, x ) = e~ iMn (x, x ) may be viewed as a diffusion 
with the diffusion coefficient proportional to the power 2 — 7 of the distance between 
the particles. When particles separate they diffuse faster and faster which results in the 
super-diffusive behavior with mean distance square growing proportionally to t 2 / 1 . On 
the other hand, on small distances the diffusion is slow and particles which get close 
spend relatively long time together. Since (/) 0x = /(xo), the time t after which 
(/). reaches, say, twice its original value behaves like O(/(x )^), i.e. it goes slower 

5 K.G. thanks E. Balkovsky, G. Falkovich and V. Lebedev for a discussion of this point 
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to zero with the diminishing separation between the initial points than for the standard 
diffusion at 7 = 2. 

We have seen that the scaling zero modes of M n correspond to conserved col- 
lective modes of the super-diffusion with the transition probabilities P n (t, x; to, Xo) . 
Existence of such conserved modes is nothing exceptional. They are present already in 
the standard diffusion. For example, 

J l%i2 - ^13] e<A ( x > x o) dx = x 2 12 - x 2 13 (4.7) 

and is time independent, although / xfx e* A (x, x ) <ix behaves like 0(t). Another ex- 
ample is 

J [ X 12 X 34 — 2(d+2) ( X 12 + X 34,)] e ( X 5 X o) <^ x = X 0,12 X 0,3& — 2(d+2) ( X 0,12 + ^0,34) • (4-8) 

Under symmetrization, the first conserved mode vanishes whereas the second one gives 
the zero mode of A whose (2 — 7)-perturbation dominates the 4-point function of the 
scalar for 7 close to 2. 

5 Some physics: short-distance behavior of fluid 
particles 

The zero mode dominance of the structure functions of the scalar is due to the appear- 
ance of such modes in the asymptotics of the Green functions M~ x (x, y) . Indeed, with 
the use of the Mellin transform, one may write (in the reduced space): 

M-^/L,y) = J± RnWL)* N-\a-lx,y) i^y)"*^*, (5.1) 

Reo-=-%- + J 

compare to Eq. ( |3.20|) . Pushing the integration contour more and more to the right and 
using Eq. (|3.18|) to control the residues of the poles, we obtain for y 7^ the asymptotic 
large L expansion: 

M-^x/^y) = j;^/,(x)^ (5-2) 

i 

with Re<7j > — 4f + | or, as we expect, with Rea^ > 0. Although it has a similar 
form to the eigen-function expansion of an operator with discrete spectrum, it has little 
to do with the spectral decomposition of M' 1 . Since M n is a positive operator in 
L 2 (R dn ) with continuous spectrum coinciding with the positive real line, the spectral 
decomposition of M~ l is a continuous integral involving the generalized eigen-functions 
of M n . The scaling zero modes /j or g; L of scaling dimensions <7j and — d n + 7 — Wi, 
respectively, are square-integrable on S' '™ -1 but are not generalized eigen-functions 
of M n (except for cr, = 0). They are rather analogous to resonances in many-body 
problems with the plane of complex a replacing that of complex energies and o with 
real part equal to — %- corresponding to real energies^. Note that due to the hermiticity 

6 this analogy is somewhat loose, since the poles in a live in the first sheet, probably only on the 
real axis 



13 



and to the overall scaling of the Green function M n 1 (x, y ) , 

M-^Ax^y) = A^ACW), (5-3) 
the expansion ( |5.2|) may be also rewritten as 

M~\L^y) = L- d ^-^ ^(x) Jfrj (5.4) 

i 

so that the zero modes gi with the scaling dimensions — d n +j— Wi of real part less than 
— ( or even < — d n +^) dominate the large distance behavior of the Green function 

of M n . Expansion ( |5.4|) may be also obtained directly from Eq. fl5.1|) by pushing the 
cr-integration contour to the left. 



It is not difficult to see directly that the functions fi appearing in expansion ( |5.2j ) 
have to describe scaling structures conserved by the flow. Indeed, 

J M-^x/^y) e-™»(x,Xo) d'x = £ IT** 'gW) J /,(x) e^ M "(x,x ) <fx (5.5) 

if we insert expansion ( |5.2| ) into the left hand side. But on the other hand, 

9 r |M- 1 (x/L,y)e- rA/ "(x,x )rf / x = d T J M" 1 ^, Ly) e^ M "(x, x ) d'x 

= -L d "^e- TA/ "(x ,Ly) (5.6) 



where we have used the scaling property (|5.3[). The (reduced space) heat kernel on the 
right hand side decays in L faster than any power for y / 0. Comparing the latter 
expression to relation (|5.5| ), we infer that d T J /i(x) e~ rMn (x, x ) tfx has to vanish and 
hence fi , a function with scaling dimension <7j , is conserved in mean by the Lagrangian 
flow. This gives another proof of the statement (4.3). 

Since the Green function is given by the time integral of the heat kernel, 

oo 

M-W) = /e-* M "(x,y)rft, (5.7) 



one may also expect to see the zero modes in the asymptotic behavior of the probabilities 
P n {t, x;0,x o ) = e~ tMn (x, x ) . Assume an asymptotic expansion 



P n (t, x/L; 0, xo) = £ L~ Pj ^(x) fy(t, x ) , (5.8) 

3 

with Repj > 0, describing asymptotics of the probabilities that the Lagrangian trajec- 
tories will come at time t close together. We could expect that functions 4>j are again 
zero modes of M n . To verify whether this is the case, consider the integral 

J P n (t, x/L; 0, y) e~ rM -{^ x ) d'x = £ L"« 4^y) j <^(x) e^ M "(x, x ) d'x. (5.9) 
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The left hand side may be rewritten as 

L d n j e -tA/„( X;y ) e - TM «(Lx,x ) d'x = / e-* A/ "(x,y) e~ L ' rM "(x, x /L) d'x 



-(t+L- r)M n(xo/L ^ y) = £ L-» 0,(x o ) ^(t + L-^r, y) (5.10) 



by changing variables x Lx and using the scaling relations (|4.5| ), the composition 
law of heat kernels and, finally, the expansion (|5.8|) . But ipj(t,y) should be smooth in 
t for t 7^ 0, oo. It then follows that 

J P„(t,x/L; 0,y) e^ M "(x,x ) d'x = ]T J df^fcy) ^(xq) . (5.11) 

j 

p=0,l,... 

The right hand side becomes independent of r only approximately if rL~ 7 <C 1 . Com- 
paring the asymptotic expansions ( |5.9| ) and ( |5.11D , we infer that the scaling functions 
<pj (of scaling dimension pj) are not necessarily preserved in mean by the Lagrangian 
flow. Instead, / 0j(x) e~ rM "(x, xo) d'x is a pure polynomial in r. Note that the poly- 
nomial still grows slower than the super- diffusive growth r Pj//7 since its order p satisfies 
the relation pj = pj> + 7P > 7P where pj> is a scaling dimension of some other </y . 
Hence functions (pj describe slow collective modes of the super-diffusion. 

It is not difficult to see how the slow modes are related to the zero modes of operators 
M n . Differentiating Eq. (|5.8|) over t we infer that 



= J2 L-p^ M n 0,(x) Vi(<,Xo). (5.12) 

j 

It follows that if the function (pj appears in the expansion then also M n <pj does. 



Since the scaling dimension of M n (pj is pj— 7 and only dimensions with real part positive 
may appear, subsequent application of M n to a (pj must produce a homogeneous zero 
mode of M n after a finite number of steps. Hence functions (pj must be organized into 
towers of descendants (pi tP based at zero modes fi = <f>i t o of M n and satisfying the 
chain of equations 

M n (pi >p = <f>i,p-x, p=l, (5.13) 

The scaling dimension of (pi tP is (cTj + jp). Since Mv +l (p itP = 0, it follows that the 
(p + l) th time derivative of 



x) e- rM "(x,x ) d'x (5.14) 



vanishes so that the above integral is a polynomial in r of degree p, in accordance 
with the previous reasoning. Note that Eq. ( |5.12|) implies that the functions ipi iP cor- 
responding to (pi iP satisfy 

i>i, P = -d t ipi, p -i = Mn^p-!. (5.15) 
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Summarizing: the asymptotics of probabilities of the Lagrangian trajectories to get 
close together is dominated by the towers of slow collective modes of the super-diffusion: 



P n (t,x/L;0,x ) = J2 £^" 7P <Mx) ^, P (t,xo). (5-16) 



p=0,l, 



Since P n (t, x; 0,x ) = P n (t, x ; 0,x), expansion ( |5.16| ) may be also rewritten as 



P n (t,x;0,x /L) = ]T L ~ ai ~ 1P M^x) <M x o) ■ (5-17) 

i 

p=0,l,... 

giving the asymptotics of the probabilities of the Lagrangian trajectories starting very 
close. The leading term on the right hand side is equal to ^o t o(t,x) — e _<Mn (x, 0) and 
it corresponds to the constant zero mode 0o,o = fo = 1. 

The above results have an important physical significance for the dynamics of the 
scalar. Recall the expression (4.4) for the time-dependence of the mean value (/) t>x 
of a function / of scaling dimension a. Inserting the relation ( 5.171) to the right hand 
side of Eq. ( |4.4| ) we obtain the asymptotic expansion of (/) t xo for large t: 

(/) t ,x = E(^" P ^W//(x)^d'x. (5.18) 

p=0,l,... 

The leading term corresponds to the constant zero mode. This term dominates for 
large t if J /(x) e~ rMn (x, 0) (fx ^ 0. However if / = <f>i >q then, as we have seen above, 
(<f>i, q ) t>Xo is a polynomial of order q in t. Consequently, / </>/ )9 (x) ^^(t, x) (fx has to 
vanish unless °"'~ cr ' + g — p is an integer between and q. Hence (f) t XQ for / equal 
to a slow mode <j>i q with a positive scaling dimension is dominated by subleading terms 
on the right hand side of Eq. ( 5.18|) . 

a_ 

To see the lower order termsQfor a generic scaling function / for which (f) tx ~ t 7 , 
it is enough to compare the mean values (/) tx for two different x . For example, 
subtracting (/) t]X for two values of x ,i gets rid of the contribution of the constant 
mode. Denote by fi yom , y ' the operator which performs the subtraction on functions 
h of x 0im : 

5 y , m ,y' 0>m h = %o,m) - %o,J- ( 5 - 19 ) 
Subtracting subsequently at two different values of x , m for m = l,...,n — 1, and 

71-1 

setting n 5 y , m ,y' 0m = S y,y'> we obtain 

m=l 

W</>v = E' (SyjM I /(x)^, P (T,x)rf'x (5.20) 

p=0,l,... 



7 K.G. thanks M. Vergassola for the discussion of this point 
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where the primed sum omits the contributions of the slow modes which do not depend on 
all variables. Hence the non-constant slow modes dominate the relative motion of groups 
of Lagrangian trajectories starting from different initial configurations. In particular, 
the relative motion is slower than the super-diffusive spread of the trajectories. This 
supports the interpretation of the slow modes as resonance-type objects in the motion 
of Lagrangian trajectories. The slow modes depending on less variables correspond to 
resonances in fewer-particle channels which drop out under the subtractions. 



6 Some mathematics: structure of the multi-body 
operators M n 

To understand analytically the origin of the asymptotic expansion ( 5.16|) , let us examine 
closer operators M n . We shall work in the reduced space TC = L 2 {H dn ). M n is a 
positive, unbounded, self-adjoint operator in Ti. Let 

(W./)(x) = e srf "/ 2 /(e s x). (6.1) 

Operators U s form a unitary version of the 1-parameter group of dilations in Ti with 
the self-adjoint operator D n of Eq. (|3.14 ) as its generator: 

U s = e isDn . (6.2) 

U s preserve the domain of M n and 

U s M n U- x = e"^M n . (6.3) 

Denote by X n the natural logarithm of operator M n : X n = In M n . X n is an un- 
bounded self-adjoint operator on Ti with the domain invariant under U s and the whole 
real line as the spectrum. The relation fl6.3| ) is equivalent to 



U s X n U; 1 = X-js, (6.4) 

i.e. to a strong form of the canonical commutation relation 

[D n ,X n ] = i 1 . (6.5) 

Since under the Mellin transform (|3.15|) D n becomes the multiplication operator by 
}(cr + 4f), X n must be unitarily equivalent to jd a by virtue of the von Neumann 
Theorem on representations of the canonical commutation relations. More exactly, there 
exists a one-parameter family U n (a), Recr = — %, of unitary operators in L 2 (S dn ^ 1 ) , 
unique modulo a right multiplication by a a-independent unitary operator, such that 

(X n fT(a, • ) = 7 Un{cr) d a U-\a) f(a, ■ ) (6.6) 
for Re a = — 4f . Equivalently, 

(M n /)>, • ) = U n (a) e^ U-\a) f(a, ■ ) (6.7) 
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or, noting that the operator d a corresponds to the multiplication by — In R n in the 
language of the original functions, 

M n = UnR^U- 1 (6.8) 

where 

(U n fT((j, •) = U n (a)f(a,-). (6.9) 

This is the promised structural result about operators M n . For the heat kernels, we 
obtain 

e -t Mn = Un(i -tn^ v -\ ( 610 ) 

What is the relation of the expression (|6.8j ) to the representation 

M n = R~^ 2 N n R-^ 2 (6.11) 

used before, with N n an operator commuting with D n l The comparison of the two 
equations gives 

N n = R2 /2 U n R-~< If' 1 R?J 2 or AT 1 = R~^ 2 U n R^U' 1 R^ /2 . (6.12) 

Let us suppose that the family of operators U n (a) has a meromorphic continuation to 
the complex plane of a with no poles in the strip 

-^<Rea<-^ + f (6.13) 

(this will prove consistent with our zero mode analysis). Then -R^ 7//2 U n R^ 2 becomes 
under the Mellin transform the operator 

e^^U n (a)e^ d " = U n (a + J) . 

The unitarity of U n for Re a = — 4f implies that 

U n (a) U n (-d n - a)* = U n (-d n - a)* U n (a) = 1 (6.14) 

and that £/ _1 (cr) also possesses a meromorphic continuation. Operator R"^ 2 U~ l R" 1 ^ 2 
becomes U~ l (o — %) under the Mellin transform and Eq. ( |6.12| ) gives rise to the relation 

N-\a-D = U n (a)U~ l (a- 7 ). (6.15) 



If N n 1 (cr — |) has a pole, see Eq. ( |3.18| ), then the simplest possibility is that either 
Un{&i — 7) is regular and then 



UJa) 



\fi){9i\UM-l) + 0{l) 



(6.16) 
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when a — ► crj or U n (<Ji) is regular and 

Unfa) = ^\fi){9i\U n {a-l) + 0(<r-<Ti) (6.17) 

with (gi\ U n (a — 7) = 0(cr — erj). In the last case, multiplying by U n (—d n — cii)* from 
the left and by U n (—d n — a + 7)* = t/ _1 (cr — 7) from the right hand side and taking 
adjoint s, we obtain 

U n (-d n -a + j) = =±=- \ gi )(fi\ Uni-dn-Wi) + 0(1), (6.18) 

i.e. relation (|6.16|) for <7j replaced by — d n + 7 — Zfj and corresponding to the twin zero 
mode Qi of /j. We expect the behavior ( |6.16| ) if ft is less singular at the origin than 
Qi and the behavior ( |6.17| ) in the opposite case (in Appendix A, this is established for 
U 2 ). Rewrite Eq. (|6TT5|) as 

Un(a + 7 p) = iV" V + 7 (p - i)) ^(cr + 7 (p - 1)) . (6.19) 

Assume that N~ 1 {ai + 7(p — J)) is regular for p = 1, 2, . . . , (i.e. that there are no 
zero modes of M n (square-integrable on S^^ 1 ) with scaling dimensions differing by 
multiplicity of 7. This should hold for generic 7. From Eq. ( |6.19| ) for p = and from 
relation (|6.16|) we infer that 

U n (a + 1 ) = -±-N-\a+l) l/<)&l ^-j) + 0(1) 

= ^ \<f>i,i)(9i\ Un{°i-l) + 0(1) (6.20) 
for a — > &i. By induction on p, it follows then that 

U n (a + 1P )) = ^ \4>i, p ){gi\ U n (<Ti-i) + 0(1) (6.21) 

where 

<P hP = iV- 1 (a, + 7 b-|)) hp-i- ( 6 - 22 ) 

Note that Eqs. ( |6.22| ) may be rewritten as the chain of relations ( |5 . 1 3| ) for the scaling 
functions cj) itP (x) = i?^ +7P 0j iP (x) where </> i)0 = fi is the zero mode of M n of scaling 
dimension <jj. By virtue of the assumption that there are no zero modes of M n of 
scaling dimension cr, + jp , the tower of descendants 4> ijP is uniquely determined^] for 
each zero mode /j. 

For / a test function vanishing near the origin, Eq. ( p.lOQ may be rewritten as 

(e-* M "/r(^x) 

00 

= JdyU^a-^JdW-^e-^'J g A CT ' y) ■ (6.23) 

R C(7 ' = _% 



^for non-generic 7 the situation may be slightly more complicated with mixing of different towers 
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After shifting the a'-integration contour infinitesimaly to the left we may perform the 
A-integral. The inverse Mellin transform of the resulting expression gives 

(e^/)(x/L) 

= ^jti L ~ a j d y Un(v,%9) / £ V{^f) (fcVlKj). (6.24) 

By moving the a-integration contour further and further to the right[], we obtain from 
Eq. ( |6.24j ) the asymptotic expansion 

(e-' M "/)(x/L) = L-^<P hP (Z) J- (6.25) 

i 

p=0,l,— 

where the sum is over scaling dimensions of zero modes of M n satisfying Reo~j > — 4f 
and where 

■ J f^~ P Y{^+ V ) (^VK^y). (6.26) 

Eq. (|6.25| ) is an integrated version of expansion (|5.16|) (at least for 7 close to 2, there 
are no scaling zero modes square-integrable on S' 1 "" 1 with — d n < Re<Tj < and, as 
mentioned before, we expect this to hold for all positive 7). 



7 Lagrangian flow 

There are some subtle points in the above discussion of the motion of fluid particles. 
When stating that, for m, n = 0, P n (t, x; 0, x ) is the joint p.d.f. of the differences 
of the endpoints of n Lagrangian trajectories, we have silently assumed that such 
trajectories, or at least their differences, make sense as random processes defining the 
Lagrangian flow on the probability space of v's. A straightforward consequence of such 
an assumption is the relation 



P n (t,x;0,x ) 



where x' = (xi, . . . ,x k ) and similarly for Xq. Eq. ( [7.1[ ) expresses the elementary prop- 
erty that the joint p.d.f. of coinciding random variables is concentrated on the diagonal. 
In particular, P n (t, x; 0, 0) should be proportional to the <i n -dimensional delta-function. 
But the heat kernels e~* Mn (x, xo) do not have this property at least for 7 close to 2 
and, expectedly, for all 7 > 0. Instead they are regular when xo — ► 0. How exactly 
P n (t, x; 0, x /L) fails to become the delta-function when L goes to infinity is described 

9 note that the poles of of the T-function do not contribute 
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by the asymptotics ( [5.17D dominated by the slow collective modes of the stochastic 



evolution of the Lagrangian trajectories. Hence, even if all joint p.d.f.'s P n of the dif- 
ferences of Lagrangian trajectories make sense as given by the m, k — heat kernels 
e -tM n ^ x ^ ^ Qijffgrgncgg f Lagrangian trajectories do not exist as random processes 
for 7 > 0. Note that for k positive we should not expect the behavior (|7.1|) since the 
Brownian motions starting from xo,« are different for different i's even if they wiggle 
around the same Lagrangian trajectory. The system behaves as if the wigglings were 
present even for k = (see more on that below). 

The 7 = and m, k = case will be analyzed in Sects. 8 and 9 and was previously 
considered in || to [TTJ, see also 0. The p.d.f. P 2 (t, x; 0, Xq) = e~ tM2 (x,xo) of the 
difference xi 2 (t) = x(t) of two Lagrangian trajectories may be easily computed in this 
case and the result is the log-normal distribution || [JTU 



P 2 (t,x;0,x ) = -X^e-^ 1 *^-^ 2 k D ^ t (x,x ) (7.2) 

where r = \x\, x — x/r and similarly for ro, xq. kt(x,xo) denotes the heat kernel on 
the unit sphere in d dimensions and it drops out in the rotationally invariant sector. 
The most important consequence of Eq. Q is that \ In f- is a Gaussian variable with 
covariance 2D/t tending to zero at large times and with mean Dd. The mean gives 
the Lyapunov exponent i.e. the rate of exponential growth in time of the distance r 
between the Lagrangian trajectories. Note that 

J r a P 2 (t } x;0 } x ) dx = r ° e D<j{ ~ d+C7)t (7.3) 

which should be contrasted with the super- diffusive behavior for 7 > described by 
Eq. (|4.6|). More generally, 



fix) P 2 (t,x;0,x ) dx = ff(e u r x) ^J=e-^ (u - tDd)2 k n ^,(x,x )dudx. (7.4) 

I \ / A.TT I ) T A 1 k 



For any test function / and for fixed t, the right hand side tends to f(0) when ro — > 0, 
in accordance with the relation (|7.1| ) and unlike for 7 > . As it is easy to see from the 
above integral (or from Eq. (|7.2|) ), the concentration of the p.d.f. P 2 (t, x;0,xo) within 
r<r] is visible if rj ^> e* Tq. For the later use, note that for a small but non-zero ro 
and for a rotationally invariant test function /, 

f(x)P 2 (t,x;0,x )dx = [ f(e tu r )^=e-^^- Dd ^ du - 0. (7.5) 



The approach of 0[|TOj was based on the observation that at 7 = the 2-point 
function of the velocity differences becomes 

( (v a (t u x x ) - v a (t u x 2 )) (v p (t2, x[) - vffo, 4)) ) = 2^5(t 12 ) 

■ [(d+ 1) 5 a/3 x 12 ■ x' 12 - x« 2 x'f 2 - x'l 2 x{ 2 \ . (7.6) 
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In particular, the first derivatives of v have space-independent correlations. In other 
words, we may set 

v(t,xi) -v(t,x 2 ) = Xtyxw or d 7 v a (t,x) = X ary (t) (7.7) 

where X a/3 (t) is a Gaussian process with values in traceless matrices with mean zero 
and the 2-point function 

(X^ a (t)X 5f3 (s)) = 2^ r <5(*-s) [{d+l)S°^P 6 - 5^5^ - 5 aS 5^} (7.8) 

obtained by differentiating twice the right hand side of ( |7.6| ) . It is easy to check directly 
that the above covariance is positive and that it is invariant under the adjoint action 
of 0(d), i.e. that X and kXk~ l have the same covariance for k e 0(d). Eqs. (|7.7f) 
are equalities between Gaussian processes. Physically, they mean that for m = and 
7 = the velocity flow acts as a uniform, volume preserving strain and rotation, as far 
as the relative motions of fluid particles are concerned. The difference of two Lagrangian 
trajectories xu(t) = x(t) should satisfy the linear (stochastic) ODE 

dx = X(t)xdt, x(0) =x (7.9) 

with a solution given formally by 

x(t) = g t ,t xo (7.10) 
where g t)to is the time-ordered exponential of an integral of independent matrices, 

— / X(s) ds ✓ _ , , N 

g tM = T e J *o , (7.11) 
of the type similar to the ones that appears in the theorems on products of independent 



equally distributed matrices |19| or in the one-dimensional Anderson localization 
The point is that g t)to may be defined as a random Markov process (a diffusion) with 
values in SL(d). It has three basic properties: 

1. gt 2 ^ and gi 2+T>il+T have the same distribution, 

2- gt 2 ,t 1 9t 1 ,t = gt 2 ,t a - e - ; 

3. g t)to is independent of g t i )t > o if (t , t) D (t' , t') = 0. 

To define such a process, it is enough to give the (transition) probability distributions 
Pt-t Q (g)dg of g t>to (dg denotes the Haar measure on SL(d)) satisfying the composition 
law: 



Pt(g)Ps(g 1 h)dg = p t+ s(h). (7.12) 

The 5'0(o?)-invariance of the Lie- algebra- valued process X imposes also the relation 

Mkgk- 1 ) = Pt (g). (7.13) 
In Sect. 8 we identify p t with the heat kernel of a certain operator on SL(d). 
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The net outcome of that analysis is that for 7 = 0, unlike for 7 > 0, the differences 
of Lagrangian trajectories Xij(t) = gt,t x o,ij are we U defined random variables. In 
particular, the knowledge of p t is all what is needed to compute the joint p.d.f.'s of 

J P n (t, x; 0, x ) /(x) d'x = J p t (g) f(gx ) dg (7.14) 

SL(d) 

for translationally invariant /. In fact, the above integrals uniquely determine pt- 



One of the consequences of the relation ( |7.14| ), closely related to the property ( |7TT 



is that, for 7 = 0, the stochastic evolution of the scalar T defined by the m, k = 

— (3 J T 2 

flow preserves the Gibbs measure formally given as ,. ® T . Indeed, the 2n-point 

1 J normalization r 

function of the scalar in this measure is 

J^ lbbs (x) = J^ pn S 8(x 12 )5(x u ) ■ ■ ■ 8{x 2n -i,2n) (7.15) 
(the odd functions vanish). But Eq. (|7.14|) implies the relation 

J P 2n (t, x; 0, xo) ^r bS (xo) rf'x = ^ GibbS (x) (7.16) 

i.e. the time invariance of the Gibbs measure correlations for 7 = 0. This should be 
contrasted with the behavior for the 7 > case where the flux of the scalar energy 
towards high wavenumbers destroys the invariance of the Gibbs measure, see ^TJ. For 



7 = 0, the invariant Gibbs measure is nevertheless unstable under perturbations, as 
follows from relation ( [7.5|) . It has also little to do with the k — ► limit of the stationary 
state of the scalar obtained in the presence of large scale forcing. The latter will be 
constructed in Sect. 9. 



The mathematics of the difference between the 7 = and 7 > cases is simple. 
Eq. (|2.8Q requires that v (t, x) be Lipschitz in x for the uniqueness of solutions^. But 
the Gaussian f-measure with 2-point function ( ^.4[ ) lives on v's which are Holder in 
x with exponent (2 — 7)/2 (modulo logarithmic corrections) but not Lipschitz, except 
for 7 = where the velocity differences become smooth, as we have seen above. Hence, 
one should not expect uniqueness of Lagrangian trajectories even if the probabilistic 
description of them may be maintained but with violation of the property (|7. 1|) . Phys- 
icallyQ, the velocity covariance should be smoothed on the dissipative scale r\ due to 
viscous effects so that it behaves as ~ Dt] -1 r 2 for r -C f], i.e. like the 7 = covariance 
with D increased to Dt]^ 1 . The Lagrangian trajectories diverge now exponentially in 
time as long as their distance is rj. Note, however, that for arbitrary small but fixed 
r o = \xo\ one never sees concentration of the p.d.f. P 2 (t, x; 0, x$) on scales smaller than 
77 if r]<e°( tri 7 - ) r , i.e. for rj sufficiently small. This explains in more physical terms 

2. 

10 recall the existence of two solutions with vanishing initial condition: x = (it)~' and x — 0, for 
1-1 

the equation x — x 2 

11 K.G. thanks G. Falkovich for a discussion of this point 
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why relation (|7.1|) fails when 77 — ^ for 7 > 0. The exponential divergence of tra- 
jectories closer than rj makes it impossible to maintain the concept of (differences of) 
individual trajectories in the inviscid limit rj — > 0. Instead, we should talk about the v- 
dependent p.d.f. P n (t,x; 0,x |t>) whose average over the velocity ensemble reproduces 
P n (t, x; 0, Xo) . It is worth noting that for positive diffusivity k, when the deterministic 
equation ( |2~8| ) should be replaced by the stochastic ODE ( |2.9| ), although the problems 
with the non-uniqueness of the solutions persist, there exists a rigorous probabilistic 
treatment^ allowing to define uniquely the transition probabilities P n (t, x; 0, x \v) for 
Holder continuous velocities |22| . Our analysis calls for an extension of such a treatment 



to the k = case. 



8 Advection by smooth velocities and harmonic 
analysis 

When 7 = and m, K — 0, the Kraichnan model becomes exactly solvable as we will 
show now. That simplifications occur in this case has been noted before, see || and j7| 
to [l3j . Our analysis is based on some observations by Shraiman and Siggia 1| . As 



was noted in [H] and in ||, for 7 = the model has extra symmetries. The operators 
M n can be expressed in terms of the quadratic Casimir operators corresponding to an 
action of the groups SL(d) and SO(d) on the correlation functions. Let us explain 
what this means. 

The group SL(d) of real matrices of determinant 1 acts on functions / on R d 
on the left by (L g f)(x) = f(g~ 1 x). The infinitesimal form of this action is given by 
4:\ t=0 L e tAf = A^Hapf where A is a traceless matrix (i.e. in the Lie algebra of SL(d)) 
and the generators H a p are 

H a p = -x a d xfj + \ 5 a px 1 d x i. (8.1) 

Similarily, on functions of n R d variables (x\, . . . ,x n ) = x, we have the (diagonal) 
action 

(L,/)(x) = /(^ 1 x) (8.2) 

with generators H a p = J2i(~~ x ?d + iS a /3x]d x i) . The quadratic Casimir of SL(d) is 
in terms of these generators 

H 2 = H a gH/3 a . (8.3) 

a,/3 

The generators of the action of the SO(d) subgroup are J a p = H a/3 — Hp a and the 
corresponding quadratic Casimir is 

J 2 = -\H J ap- (8-4) 



12 we thank G. Eyink for pointing this out to us 
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The observation of Shraiman and Siggia was that when 7, m = in the velocity covari- 
ance d af3 (fD|, the operator M n = J2i<j d a/3 (xi - : Cj)d x <*d p becomes 

M n = ^[(d+^J 2 - dH 2 }. (8.5) 



By definition of the action (|8.2| ), the same formula holds also when we express M n in 
terms of the n — 1 difference variables. In particular for M2, the Casimirs H 2 and 
J 2 correspond to the action of SL(d) and SO(d) on functions f(x) of the difference 

variable x = x^. In this case, we may diagonalize H 2 by the Mellin transform f(x) — > 

^ 00 

f(cr,x) = J r -0 " -1 f(rx) dr : 


(H 2 fr(a,x) = d -^a{a + d) f(a,x). (8.6) 

with Rea = — |. It follows, in particular, that the spectrum H 2 acting in L 2 (R d ) is 
] — 00, anc i that of M 2 is [^j-, oo[. Denoting e~ tj2 = k t , we obtain 

/(x) e-* M2 (x,x ) rfa; 

00 

Reo-=-f 

00 

= JfJdxf(rx)^ t e-^^~ tD ^k D ^ t (x,x ) (8.7) 



where we have performed the Gaussian integral over a. The result (|7.2j ) readily follows. 
Taking f(x) = C^(r) = C(r/L) with C the rotationally invariant forcing covariance, 
we obtain by integrating over t the expression for the 2-point function of T at 7 = : 

T %L {x) = (M 2 - l C L )(x) = ^ (£c(p/L) d -Z + r~ d f o C{p/L)p d ~Hp^ . (8.8) 

Clearly ^^(x) is smooth for 1 ^ and 

T %L {x) = In (r/L) for small r, (8.9) 

^2,l(x) ~ (r/L)~ d for large r. (8.10) 



In order to solve Eqs. (|2.23| ) for the higher-point functions of T we need a rep- 
resentation for the Green function M" 1 . This is obtained by relating H 2 and J 2 
to the Casimirs 7i 2 and J 2 of the left action of SL(d) and of SO(d) on func- 
tions F on SL(d), given by (C g F)(h) = i^g" 1 ^), or in the infinitesimal form by 
■^\ t=0 C e tAF = AP a l-L a pF . Note that T~i a p are skew-adjoint in the regular representa- 
tion and that 

(d + i)j 2 - du 2 = -^£(^-^) 2 - - A Y,^ + n Pa f (8.11) 

a,/3 a,/3 
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is a positive elliptic operator in L 2 (dg). In particular, it has the heat kernel 

K t (g,h) = e-^ d+l ^ 2 - dH \g,h) (8.12) 

satisfying JJC t (g,h)dh = 1 and JC t (g,h) = }C t (kg,kh) = K t {gg' , hg') for k G SO(d) 
and g' G SL(d) . Assign to a translationally invariant function /(x) and to x a function 
F*.{g) = /(px) on SL(d). The linear map f *—>■ F intertwines the two actions of SL(d): 

(L g f)(hx) = (C a F x )(h). (8.13) 

It follows that 

Je- tM "(gx,y)f(y)d'y = J JC t (g, h) f{hx) dh. (8.14) 

Comparing the above relation to Eq. (|7.14j ) we conclude that 

Mhg- 1 ) = K t (g,h). (8.15) 

Clearly the basic properties (|7.12|) and (|7.13|) of the p.d.f. p t follow. In other words, 
we may identify g t _ to as the diffusion process on group SL(d) with the generator equal 
to ^[{d+l)J 2 -dH 2 }. 

Integrating the relation ( 8.14j) over t we infer that 

j M-\g^y)f(y)d'y = J Q{g,h) f(hx) dh (8.16) 

where Q is the integral kernel of (a?r[(<i+ ^)J 2 — dH 2 })^ 1 . Applying iteratively 
identity Q8.16|) to Eq. (p.23|) we end up with the expression 

JF 2 „(x) =^F 2n (u p ) (8.17) 

v 

where the sum runs through all ordered pairings p = ({h,ji}, • • • , {i n ,jn}) of {1 . . . 2n}, 
Up — (xi 1 j 1 , . . . , X{ n j n ) and 

„ n „ n 

F 2n (u) = / YlG{gir-i,9i)C{giUi)dgi = / X\Q{gi-i,gi)C(g i u i ) dg t (8.18) 
J i=i J i=i 

where go = e and Q(g, h) = J SO(d) G(g, kh)dk. The last equality follows by substituting 
9i = hg'i, 92 = k x k 2 g' 2 , g n = k x --- k n g' n , and using Q(kg 1 , kg 2 ) = Q{gi, g 2 ) and C(kx) = 
C(x). 

The final reduction consists of identifying Q with the Green function of the Laplace- 
Beltrami operator A on the homogeneous space Hd = SL(d)/ SO(d) . By definition, 
A coincides with the Casimir ^Ti. 2 if we view functions on Hd as functions on SL(d) 
right-invariant under the action of SO(d). Assign to a function / on Hd the function 
9 - f{9) = fig- 1 )- Clearly f(kg) = f(g) and (CJ)^ 1 ) = f(hg) = (ll g f)(h), i.e. 
the map />—>•/ intertwines the action of SL(d) on the functions on Hd with the right 
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regular action of SL(d). Since the quadratic Casimirs of SL(d) in the left-regular and 
in the right-regular representations coincide and J 2 vanishes in the action on /, we 
infer that 

- 2d(Af)(g- 1 ) = -d(H 2 f)(g) = {[{d+l)J 2 - dH 2 ]f)(g) (8.19) 

and that 

J G(g-\ h) f(h) dh = J g(g, h) f(h) dh = J g(g, h) f(h) dh (8.20) 

where the function G(g, h) on x represents the kernel of the operator (— D'A) -1 
where D' = . Thus Q{jg,h) = G(g~\ /i" 1 ) and Eq. (|A8|) becomes 



/n 
UG(g t -u9 t )C(g^u t )d gi (8.21) 
i=i 

Every matrix g G SL(d) can be uniquely represented as a product (the so called 
Iwasawa decomposition) g = nak where k e SO(d), n is upper triangular with 1 on 
the diagonal and a is diagonal with positive entries. Thus one may parametrize the 

cosets gSO(d) by na. For d = 2 we may write a = diag(?/^, y~2 ) , y > 0, n = (J \ ), 

a; G R. The Haar measure <ig becomes dg = y~ 2 dxdydk. The homogeneous space 
H 2 may be identified with the upper half-plane H = {z = x + iy G C\y > 0}. 

The action of SL(2) on H is given by the Mobius transformations (" b d )z = 

Since ki = i, the identification maps the coset gSO(2) to gi = nai. We shall denote 

i _ i 

na = g(z) = (^ v The SL(2)-invariant measure on H is du(z) = y~ 2 dxdy and 
the Laplace-Beltrami operator becomes 

A = y 2 (d 2 y + d 2 x ). (8.22) 

The Green function G is given by the explicit expression: 

G(2 , ,) = ^_ ln p^l+pjlL . (8 .23) 
Eq. ( |8.21| ) may now be rewritten as 

^2„(u) = \^G{z i ^z l )C{g{z l )- x u>)dv{z i ) (8.24) 
with zq = i. In Appendix B we study the integrals (8.24) in more detail. In particular 



we show that the leading singularities at coinciding points of the correlation functions of 
T are given by a Gaussian expression, a sum of products of 2-point functions, confirming 
the analysis of [0 . 
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For the dimension d > 2 one can proceed analogously. In the Iwasawa decompo- 
sition we parametrize n by the off-diagonal entries, x a , a = 1, and write 
a = e^, 4> = diag(0i, . . . , 0^) with J2i(fii — 0. The Haar measure becomes in these 
variables 

dg = eEi<i^-^ 11 d<Pi J] dx a dk . (8.25) 

1=1 a 

G is (proportional to) the Green function of the Laplace-Beltrami operator A on 
SL(d)/SO(d). Explicitly for d = 3 write = Jadiag(l, -1, 0) + i/Mag(l, 1, -2). 
Then 

A = e 2a d 2 Xi + e a+(3 d 2 X2 + e p - a {d X3 + Xl d X2 ) 2 + d 2 a + 3d 2 . (8.26) 

and dg = e~ a ~ ,3 dadf3dxidx2dx^dk. There does not seem to exist a very explicit 
expression for G in d > 2. However, the singular behavior of jF 2n can be extracted 
again, see Appendix B. 



Let us end this section by deriving the formula for the Lyapunov exponents of the 
Lagrangian trajectories, previously found in [Tl[ by path-integral techniques. In J| it 
was observed that M n may be also expressed using the quadratic Casimir of the action 
of SL(n — 1) with the generators 

G« = -x? n d x? + (8.27) 

for 1 < i,j ' < n— 1. This action corresponds to the natural action of SL(n — 1) on the 
i-index of Xi n = Xi — x n . Denoting by G 2 the quadratic Casimir J2 i:j GijGji and by A 
the generator of dilations xf d x <*, one obtains 

M n = ^[(d + l)J 2 -dG 2 - d -^A(A + d n )}. (8.28) 

Let p denote the volume spanned by vectors X{ n , i — 1, . . . , n — 1, describing the time 
t differences of the Lagrangian trajectories starting at time zero from points x : 



p = ^det^Xin - x jn ). (8.29) 
We would like to find the p.d.f. of p. Note that for a function f(p), 

M ^f(P) = ~WWS^^ A + d ^ M = - ( - 1)( l7 +1)P P d p {pd p + d) f(p). (8.30) 

This follows from Eq. (|8.28 ) since p is SL(n — 1)- and S , 0(c?)-invariant. Hence M n 
preserves the space of functions f(p). Also 

J\f(p)\ 2 ]Jdx in = const. / imf^dp (8.31) 

4 
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where const. = / 6(p — 1) II, dx in . Hence M n in the action on f(p) is diagonalized by 

^ oo 

the Mellin transform f(p) — ► /(cr) = / p~ a ~ 1 f(p) dp (unitary for Recr = — ~): 



o 



(M n /)» = _ (n-D(^+i)^ (T(a + rf)/(a) s _ L) ( n ) (T ( (7 + d ) / ( (T ). (8.32) 



d-1 



As in Eq. ( |8.7| ), we obtain 

oo 

/ /(p)P n (t,x;0,XoKx = J%f(p) -^=e-^ {ln ^~ tD{n)d)2 . (8.33) 

o * 

Hence | In p is a Gaussian variable with covariance 2D(n)/t tending to zero at large 
times and with mean D(n)d which, by definition, is the sum of the (n — 1) largest 
Lyapunov exponents describing the effective separation of n Lagrangian trajectories. 
We infer that the n th Lyapunov exponent is 

A n = (D(n+1) - D(n))d = \{d - 2n + 1)D' (8.34) 

with d exponents equally spaced and symmetric with respect to the origin, confirming 
the result of HIT . 



9 Quadrature of the 7 = case 

Let us explicitly construct the stationary state of the passive scalar advected by smooth 
Gaussian velocity with 2-point function ( |7.6| ). Relations ( |B.17| ) and ( |8.21| ) allow to write 
a compact expression for the generating function of the 7 = theory: 

00 . 

*(X) = (e*/^*)*) = ESE / FMUxix^dx,. (9.1) 

n=0 V J i 

Noting that all pairings give the same contribution to the Xi integral and that there 
are of them we get 

00 « n 

$( X ) = £(-!)" / n^ft)^)^ (9-2) 

where 

V x {g) = \ Jc(g- 1 (x~y))x(x) X (y)dxdy (9.3) 

is a non-negative function on H d bounded by V x (e) = ?C(Q) (J x) 2 ■ Eq. fl9.2|) may be 
rewritten in the operator language as (D f = Wj) 

00 . 

HX) = EC- 1 )" /[(^^) n l](e). (9.4) 

n=0 J 
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The sum on the right hand side involves the Neuman series for the operator (—D'A + 
V^) -1 , i.e. for the Laplacian on H d perturbed by a potential. Resumming the series we 
obtain 

= 1 - [(-D'A + V x )- l V x }(e) (9.5) 

which is an explicit expression for the characteristic functional of the stationary state 
of the 7 = Kraichnan model. 

Let us see that the right hand side of Eq. ( |9.5|) makes sense. Using the Feynman-Kac 
formula expressing the perturbed heat kernel as an expectation E g (-) with respect to 
the Brownian motion on with transition amplitudes e~ D A (g, h) , starting at time 
zero at g: 

e -t(-n>A +Vx ) igjh) = Eg (e-Io v ^ ds 5 h (h(t))y (9.6) 
we infer the bounds 

0<e-«- D ' A+v x\g,h)<e tD '*(g,h), (9.7) 
< (—D'A + V x )-\g, h) < G(g, h) . (9.8) 

Since 

[(-D>A + V x )-%](e) = j(-D'A + V x y\e,h)V x (h)dh, (9.9) 
it follows that the latter integral is bounded by the smeared 2-point function 

J G(e,h)V x (h) dh = \J -^(ziz) x{ x i)x{ x 2) dxxdx 2 (9.10) 
which is finite for test functions x e -g- from the Schwartz space iS(R d ), see Eq. (|8.8| ). 

$ defines a continuous positive-definite functional on iS(R d ). The continuity of 
$(x) w.r.t. x £ <S(R d ) is easy: it follows by the Dominated Convergence Theorem 
from the Feynman-Kac representation of the perturbed Green function: 

(—D'A + V x )-\g,h) = J~dt E g (e-fo v - Ws))ds 5 h (h(t))^ . (9.11) 

The positive definiteness: 

J2X r T s ^(Xr-X s ) > 0, (9.12) 

r,s 

is a little bit more complicated. Let us sketch its proof. Define first the positive definite 
characteristic functional 

$ t ( X ) = ( e i/T(t,*)x(*)«fc) ( 913 ) 
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of the time t (quasi-Lagrangian) state of the scalar where T(t, x) = J * f(s, g^x) ds is 
a functional of the forcing / and of g ttS . The above expression for T(t, x) is obtained 
for the initial condition vanishing at t = in Eq. (|2.4j ). The expectation in ( 9.13| ) 



is w.r.t. the Gaussian measure of the forcing and w.r.t. the measure of the diffusion 
process g tyS - It is easy to see that J T(t, x) dx is square-integrable with respect to 
these measures. Performing the integration with respect to /, we obtain 

The remaining expectation over g t a is easy to calculate by expanding the exponential 
(the resulting series of expectations converges absolutely for finite t). The result is 

MX) = I" f[e- s ^ D '^V x ]{e)ds. (9.15) 

Using the bound (|9.7|), it is easy to see that $ 4 (x) converge to $(x) when t — > oo. 
Hence the positive definiteness of Note that Eq. (|9.15|) may be rewritten by inte- 
gration by parts and the Feynman-Kac formula ( |9.6| ) as 



$ t ( x ) = [e-'(- D ' A+ ^)l](e) = E e (e-Io v ^ h ^ d A (9.16) 

which follows also directly from Eq. ( |9.14| ) if we notice that the diffusion process s i— > 
gt,t-s on SL(d) projects to the Brownian motion on Hd- The resulting alternative 
expressions for 

= lim |e-'(- D ' A+ ^(e,/i)^ = E e L~Io°° v ^) ds ^j . (9.17) 

relate $(x) to the long time behavior of the diffusion on the homogeneous space iJ^ in 
the presence of a positive potential V x or to the low-energy properties of the Schrodinger 
operator —D'A + V x . They imply that < $(x) < 1- Expressions (|9.17| ) may also 
be obtained directly in the Martin-Siggia-Rose (MSR) ]23| formal functional integral 
approach. 

By Minlos Theorem, the normalized ($(0) = 1), continuous, positive-definite func- 
tional $ on (S(R d ) given by Eqs. ( |9.5|) or ( |9.17|) defines a unique probability measure 
dfi on S'(R d ) s.t. 

$( x ) = [ e'l^^^d^T). (9.18) 



dfi is the stationary state of the Kraichnan model for 7 = alluded to in Sect. 7. 
It is quite different from the Gibbs measure and quite non-Gaussian and is, indeed, 
supported by distributional configurations of the scalar since the correlation functions 

JT 2n ( x ) = Jt{x x ) ■ ■ ■ T(x 2n ) dn{T) (9.19) 

diverge logarithmically at coinciding points. The measure dfi contains all the joint 
p.d.f.'s of smeared scalar values J T(x) x( x ) dx. In particular, the function p \— > $(px) > 
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is the Fourier transform of the p.d.f. p x (0) of / T(x) x( x ) dx whose behavior was studied 
in [0, see also 0@|ig . 



$(px) 1S a pointwise limit of the finite-time functions $t(px) which are entire in p. 
For Rep 2 > -b 2 , b > 0, 

|$ t (px)| < $ t (±%) = E e (e b2 fo v ^ h ^ ds ') = J e^ D ' A+b2v -\e,h)dh 

= 1 + b 2 JdsJ e s ( D ' A+b2v *\e, h) V x {h) dh. (9.20) 

The Schrodinger operator —D'A — b 2 V x with a negative potential may develop bound 
states. The right hand side of the inequality ( |9.20| ) grows with t since the expression 
under the integrals is positive. If e& = inf {spec(— D'A — b 2 V x )} < then the growth is 
unbounded since e s( - D ' A+b2Vx \e, h) ~ e~ SSb for large s. On the other hand, for > 
the right hand side of ( |9.20| ) would be bounded uniformly in t if V x were of compact 
support on Hd- V x , however, does not have a compact support as a function on Ha 
even if C and \ do (if they do not vanish identically). It is, nevertheless, easy to see 
from the definition (|9.3|) that V x vanishes at infinity of H^, i.e. that it gets arbitrarily 
small outside sufficiently big compact subsets of . This is enough to assure a uniform 
bound for the right hand side of ( |9.20D as may be seen by the following argument which 
separates the behavior at infinity of Hd from that in the interior. Write V x — V' + V£ 
where < V' < V x and V' has a compact support. By the Holder inequality, 

If we choose e small so that for b' = (1 + e) 7 6 the relation ey > still holds then the 
first expectation on the right hand side of inequality ( |9.21| ) is bounded uniformly in t (e& 
increases with decrease of V x ). Choose the support of so that ^b 2 V^ < v o < 
where vq is a constant. Then 

E e (e^ b2 Jo W))*^ < 1 + ^T b2 f o ds j e~ s(D ' A ~ Vo) (e,h) V x {h) dh 

= 1 + ^b 2 (ds [ e- s ^~ V0 \x 12 - y) C(y)x(x 1 ) X (x 2 ) d Xl dx 2 dy (9.22) 



and the last expression is bounded uniformly in t as may be easily seen from Eq. ( |8.7| ). 

We infer that $(px) > as a limit of uniformly bounded analytic functions, is analytic 
in p for Rep 2 > 6q but has a singularity at p = ±z6 where b is the positive number 
s.t. = 0. The Cauchy bounds imply now that 

I I < c -^Sw^ ( 9 - 23 ) 

in any strip |Imp| < bo — e for e > 0. Note also that lim ${px) — by virtue of 

Ip|-*oo 

Eq. ( |9~T7| ). Since 

p x {9) = ^Je- ie ^(p X )dp, (9.24) 
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it is easy to show integrating by parts few times and moving the p-integration contour 
to Imp = ±(6 — e) that p x {0) is smooth except, possibly, at 9 = and that 

p x (6) < const. (e) e ~ {b °- t)W (9.25) 

for \6\ > 0(1) and any positive e. Clearly, the same inequality fails for negative e 
since it would imply analyticity of $(px) at p = ±ibo- I* 1 short: the p.d.f. p x (0) 
of / T{x) x( x ) dx has an exponential decay for large 16*1 with the rate bo equal to the 
value of b at which the ground state of —D'A — b 2 V x crosses zero energy. Note that the 
rate bo, as related to a bound state energy is not, in general, a semi-classical quantity. 

For rotationally invariant \ which, for simplicity, we shall normalize so that J \ = 1 , 
our operators on H d reduce to the ones on the double coset space SO(d)\SL(d) / SO(d) . 
This space may be identified with the Cartan algebra of SL(d) divided by the action 
of the Weyl group and may be parametrized by the diagonal matrices diag(0i, . . . , 4>d) 
with entries <p% < • - • < (j>d and s.t. Yli<Pi — 0. In this parametrization, the Schrodinger 



operator —D'A + b V x becomes the Calogero-Sutherland Hamiltonian |24j [2l| with a 
potential: 

^f-jE^-Eiss^ + ^l-^W- (9 - 26) 

acting in L 2 { T7 dcpi) . The constant d<yd 2 ^ , equal to the half length squared of the Weyl 



vector of SL(d), is the infimum of the spectrum of —A [25] so that > D d ^ 4 — — — 

b 2 V x (0). Note that, for d > 2, d ^ 2 ~ 1 ^ is higher than the infimum of the spectrum of 
— \H 2 acting in L 2 (R d ) since, as pointed out in the remark after Eq. ( |8.6j ), the latter 

is equal to This discrepancy is due to the appearence of different irreducible 

representations in the decomposition of the actions of SL(d) in L 2 (Hd) and in L 2 (R d ) 
for d > 2 and it will play an important role below. When the forcing covariance C(r) is 
essentially constant for r < L and is falling off to zero for r ^> L (e.g. for C(r) replaced 

by Cl{t) = C(r/L)) then the potential —b 2 V x approaches for L — > oo a constant equal 
to -\b 2 C{Q) so that for large L we obtain hb 2 C(0) = D ' d %-v or 



io-d^. (9.27) 

Note however that although bo stabilizes when L — > oo, the right hand side of Eq. 
( ggg ) tends to e ? tb2c{0) and blows up with t. 

The exponential decay of the scalar p.d.f. for 7 = in the isotropic two-dimensional 
situation was first found in 0, see also for a discussion of the non-isotropic case. 
The calculation of was extended to higher dimensions in |l(J. Both calculations 



were reinterpreted in |T3 within the semiclassical approach. Our rigorous result about 



the decay rate bo of p x (0) disagrees for d > 2 with the result of Jl0| and with the 



instanton calculation of [|T^]. These papers obtain the value b' = d J for the decay 
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rate which is smaller than bo for d > 2. The point is that in || and [|K|] the function 
V x (g) of Eq. (|9~3| ) was replaced by V x (g) = ^C^g^x) with fixed x ^ 0. This simplifies 
the calculation of the expressions of Eq. ( |9.17|) since only the distribution of g tjS x for 
one x is needed. They become 

$' = E e (e-ST V 'M'»*') = lim / e- t{Ah + ^ c \x,y)dy (9.28) 

and lead to the quantum mechanical problem analyzed in [P| ||10|| . Upon the replacement 
of V x by p 2 V x one obtains a function ^(p) whose first singularity off the real axis is 
at p — ±ib with b s.t. the ground state of M 2 — \b 2 C crosses zero energy. Since the 
spectrum of M 2 starts from see the remark after Eq. (|8.6| ), we indeed obtain, for 



4 

C = Cl and large L , the exponential decay rate b' for the Fourier transform of <&' x (p) . 
The technical reason for the discrepancy with our exact calculation is that V x (g) , unlike 
its smeared version V x (g) , does not vanish at the infinity of and leads to a more 
singular behavior of the right hand side of Eq. (|9.20| ). Another way to see it0 is that <&' x 
is given by a version of Eq. ( |9.1| ) with / T7 x( x i) dx% omitted and with F 2n (up) replaced 
by the partition-independent contribution F 2n (x, . . . , x) corresponding the collinear 
configuration u p giving the most singular behavior when u p — > (see Appendix B). The 
smearing in Eq. (|9.1| ) makes this behavior more regular. Our result persists, however, 
also if we replace V x (g) with V^(g) — j J C(g~ 1 x)ip(x)dx, if ip and C are non-negative 
function from <S(R d ), since V^(g) still vanishes at infinity. In particular, ip may vanish 
around the origin which shows that it is the smearing of collinearity, not the inclusion 
of coinciding points, which is responsible for the discrepancy between 6 an d b' Q . The 
lesson is that the correlation of (non-collinear pairs of) Lagrangian trajectories renders 
the smeared scalar less intermittent in more than two dimensions and should not be 
neglected. 

It is easy to see that $(p%) decays exponentially for large real p. Denote by r the 
first exit time of the Brownian motion on from a fixed neighborhood of e. The 
probability of a given value of r is bounded by e~ const -/ r . Since V = f£°V x (h(s)) ds > 
const. r, the conditional expectation E e (e~ pV \r) is bounded by e ~ const p r . Hence 
the exponential decay of E e (e~ p v ) < j^° e ~ const -(p t+i/t) ^ r _ more exac t description 
of the decay follows from the path-integral integral representation of the expectation 
(|9.17|) . The latter implies that the large p behavior of $(px) f° r rea l Pi unlike the 
large 6 behavior of p x (&), is semi-classical: 

HPX) ~ e- |p|s{9( - )} (9.29) 

where [0, oo] > g(s) describes a trajectory (instanton) in minimizing the action 

S(K-)) = J^\h(s)\ 2 + V x (h(*))) ds ( 9 - 30 ) 
for fixed initial value h(0) = e (with | ■ | 2 standing for the SX(<i)-invariant metric on 



Hd). This is the same instanton as in the field theoretic MSR analysis of ]T3|]. For 
13 we thank M. Chertkov for suggesting this interpretation 
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rotationally invariant \i the problem reduces to the one on SO(d)\SL(d)/ SO(d) with 
the action 



and with the initial value 



f(iyT,Hs) 2 + 1^(0,00)) ds 



0. In d = 2 the minimal value of S is 



(9.31) 



i 



\JvMd<t> > 



(9.32) 



where 6 = do. — 



proximately constant up to 



For C(r) approximately constant up to r = L, V x (<f>) is ap- 



lnL and then it decays to zero like ~ e 



Con- 



sequently, the exponential decay rate of $(px) is approximately InL J^M- for large 



L, in agreement with and |L3[ 



For d > 2 and large L the minimum of S is at- 
tained on the trajectory which in the region of constant potential goes in the direction 

1, d — 1) and the value of the action is again 



i 



d(d-l) 



InL 



up to lower 



2D 



order terms, as pointed out in |10[ and |p~3 
that p x (0) is smooth also at zero. 



The exponential decay of 3>(p%) implies 



10 Conclusions 

In this paper we have analyzed the stochastic dynamics of Lagrangian trajectories for 
Gaussian, time-decorrelated random velocity fields considered in the Kraichnan model 
of passive advection. We found that the dynamics is characterized by two related phe- 
nomena. First, the Lagrangian trajectories loose in the limit of high Reynolds numbers 
the deterministic sense for a fixed velocity realization due to their sensitive dependence 
on initial conditions. Second, their relative stochastic dynamics is dominated by slow 
resonance-type modes. The slow modes determine the average characteristics of the 
spread of Lagrangian trajectories responsible for the loss of their deterministic charac- 
ter. Both phenomena were essentially due to non-smoothness of the typical velocities 
signaled by fractional Holder exponents in their spatial dependence. Since the turbu- 
lent velocities are non-smooth in the limit of high Reynolds numbers, we expect the 
two phenomena to persist for more realistic velocity ensembles and to continue to be 
responsible for the anomalous scaling. For the spatially smooth velocities, we calcu- 
lated the Lyapunov exponents describing the sensitive dependence of the Lagrangian 
trajectories on initial conditions for distances smaller than the viscous scale. Using 
harmonic analysis on the symmetric spaces SL(d)/SO(d) we also obtained in this case 
an explicit form of the characteristic functional of the stationary state of the passive 
scalar and exhibited an exponential decay of the p.d.f.'s of smeared values of the scalar 
relating the decay rate to the properties of the ground state of the Calogero-Sutherland 
Schrodinger operator with a potential. 



Appendix A 
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We shall make explicit the structural result of Sect. 6 for the heat kernel e tM2 (x,xo) . 
In the angular momentum I — 0, 1, . . . sector, 

M 2 = M 2 (0 = --^d r r d+1 ^d r + W£pLl(d-2 + l)r-~< (A.l) 

which is a positive operator in L 2 (]0, oo[, r d ~ 1 dr). The generalized eigen-f unction of 
M 2 (7) corresponding to eigenvalue E > involves the Bessel function 

Mr) = r^U^r 1 ) (A.2) 

where 



vi = \yj{d - 7 )2 + A d -^l{d - 2 + /) . (A.3) 
The spectral decomposition of M 2 (/) has the form 

M 2 (l) = j E\ip E ){ip E \dv{E). (A.4) 

Since 

{U.<p E ){r) = e srf/ V E (e s r) = b\^ e {t), (A.5) 

we infer that 

U s M 2 (l)U- 1 = J E\<p E )(<p E \dv(e^ s E). (A.6) 



Since, on the other hand, W S M 2 W S 1 = e 7S M 2 , see Eq. (|6.3|) , it follows that 

du{e ls E) = e~< s dv{E), (A.7) 
i.e. that dv(E) = cdE for some positive constant c. Hence 

c J \<PE)(<p E \dE = I (A.8) 

and for 



f{E) = V~c / VE {r)f{r)r d - l dr, (A.9) 

•'0 

we obtain J o °° \ f(E)\ 2 dE = J o °° (/(r)^ 6 '" 1 dr. Substituting E = e "", we shall define 

7 

(Vi/)(«) = ^^"/(e 7 "). (A.10) 



Vi : L 2 (]0, oc[, r d x dr) — > L 2 (R, rfw) is a unitary operator. Besides, 

(ViM 2 (/)/)(tt) = e*^(Vi/)(«). (All) 
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Let V 2 : L 2 (]0, oo[, r d 1 dr) — ► L 2 (R, dw) be another unitary operator defined by 

— in 

(V 2 f)(u) = e 5 /(O. (A.12) 

Note that 

(ViUJ)(u) = (Vif)(u-s), i = l,2, (A.13) 
so that ?7 2 — 1 /? f 1 V2 commutes with U s . Besides, 

M 2 (0 = U 2 r-^\ (A.14) 
as follows from Eq. ( |A.11[ ). This is the relation flt>.8| ) for n = 2. Since, explicitly, 

(Vxf)(u) = VTcJ U^B^ U>) ) (V 2 /)K) du' (A.15) 

and the Mellin transform is the composition of V 2 and the Fourier transform, we obtain 
for Re a = I 

d+2a 



■J- 1 Ux)dx = ^(iVD) T{ l {1 + Ui _^ )y (A-16) 
The unitarity implies that 'jcD = 1 so that, finally, 



The right hand side has a meromorphic continuation to the complex plane of o with 
poles at 

<ri* = + \n +1P (A.18) 

for p = 0, 1, . . . Since the true (more regular at the origin) zero mode of M 2 (/) occurs 
at scaling dimension ai >0 , this is exactly the analytic structure predicted for U n (a). 
The function r ai,p (multiplied by an angular term) represents a slow 2-point mode in 
the angular momentum I sector. 



Appendix B 

Let us briefly consider the convergence properties of the integrals ( |3.24| ). Let 
ki e SO{2) be rotation matrices s.t. Ui = fcj(rj,0) where = \m\. We have g[zi)~ x ki = 



37 



(fcj l g(zi)) 1 = k' i g(k i 1 z i ) 1 for some k\ G 5*0(2). Denoting C(k(r, 0)) = C(r), observ- 
ing that ^(zj) -1 ^, 0)| = 2 and using the SL{2) invariance of dv(zi) , we obtain 

F 2n (u) = J G(i,zi)G(K 1 z 1 ,z 2 ) . . . G(K n ~iZ n -i, z n ) JlCinyi 2 ) dv(zi) (B.l) 

i 

where = k~+ x ki is the rotation by the angle between Ui + i and Ui. We shall study 
the behavior of F 2n as rj tend to 0. The following is a useful relation: 

J G(z, z') dx' = -L - y ) + - y')) = G (y, 2/') • (B.2) 

For the 4-point function, noting that Im(«i2i) = 712/1 where 

71 = [(a;isin^ + cos^) 2 + yj[siii 2 0]~ 1 , (B.3) 
•d being the angle between u 2 and izi, we obtain 

/_i _ 1 

G(i,z 1 )G (jiy 1 ,y 2 )C{r 1 y 1 *)C(r 2 y 2 2 ) y 2 2 dy 2 dv{z x ) . (B.4) 

Consider first the case •& = i.e. 71 = 1. Then 

F 4 («i, u 2 ) = / (0(2/1 - 1) + - yi))(yie(y 2 - y x ) + y 2 %i - y 2 )) 

_ 1 _ 1 

■ C(ri2/! 2 )C{r 2 y 2 2 )y^ 2 dy l y^ 2 dy 2 . (B.5) 

Since C = Cl has rapid decay at infinity, the integrals are effectively cut to yi > (ri/L) 2 
and produce logarithms of (r^/L) as these ratios tend to zero. The most singular 
contribution is from y\9(l — yi)y 2 Q{y\ — y 2 ) term which yields 41n(r 1 /L) hi(r 2 /L) — 
2(ln(ri/L)) 2 if n > r 2 and 2(ln(r 2 /L)) 2 if r 2 > r 1 . Thus 

-^4,l(«i,m 2 ) + F A (u 2 ,Ui) = {^-) 2 ln(ri/L) ln(r 2 /L) + less singular (B.6) 

for ■d = 0. 

For 7^ 0, the y 2 integral yields 

JW«i,«2) = £ / ln !f + - |)2 ^(^i^ 1 ) [(MtO + Hvir 2 2 L 2 )) 

■ #(im - r 2 2 /L 2 ) + B] yfdy l dx 1 (B.7) 

where B is bounded. For the ln(7i) term the only singularity is at y\ small and this 
term is bounded by const. | m(ri/L)| . The rest has same leading singularity as the 
•& = calculation. Thus recognizing in ([B~6D the 2-point singularities (|8.9|) , we infer 



that the leading singularity of the 4-point function is Gaussian, the sum of products of 
2-point functions. 

The analysis of the general correlation is similar though tedious. When all the 
points are on the same line i.e. all the angles are zero, we can do all the Xi integrals 
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by ( |B.2|) . Most singular contribution is the one where all Cro(z/i-i> are replaced by 
Ab~yi8(yi-i ~ Vi) ■ Summing over the permutations of the yields 

EW=Ili [ C L {r t y^)y-My + . . . = U~^ Hn/ X) + ••• (B.8) 

where ... is less singular. Non-zero angles give again subleading contributions. The 
subsequent sum over unordered pairings gives the Gaussian expression for the leading 
short-distance singularity of J r 2n,i( x ) m terms of the singular contributions to the 2- 
point function, in agreement with the observations of || and ||12|| . 

Finally, for d > 2 , to extract the leading singularity some bounds for G are needed. 
Let us here note only that for d = 3 if all the «j are identity, then the x integrals can 
again be done and the result is that G gets replaced by Go where 

G = const, (-yldl - SyldD' 1 (B.9) 

on L 2 ((yiy 2 )~ 2 dyidy 2 ) (we have put y\ = e a , y 2 = e 13 in ( 8.26|) ). The behavior of Go 
near yi = is calculable and the leading singularity can again be shown to be given by 
products of 2-point functions. Using the Mellin transform we may write 

(GoFM, y ' 2 ) = const./ ^e-/^(^)%^ (1OS t )2 -^ (log f )2 F( yi ,y 2 ). (B.10) 

Cij-y^^) in ( |B.1| ) is replaced by C(ry 1 2 y 2 6 ). Hence, let p = y{y 2 and let F(yi,y 2 ) = 
f(p). Then (G F)(y' 1 ,y' 2 ) = f(p') with ' " 

f(p') = const, f e"* J dudv /(e^V) e u+3v e"^"^ (B.ll) 

which after performing the u — v and the t integrals becomes / g(p', p) f(p) dp with 

g{p\p) = const. [\0{pt -p) + ^6(p- p')] . (B.12) 

We may then write F 2n in the form 

F 2n {u) = I g(l,pi)...g(ft n -i,p n )Yl c ( r iPi 1 ) d Pi (B.13) 

i=l 

and the analysis of the singularities goes on as in the d = 2 case. 
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